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ABSTRACT 


V 


It  j s  assumed  that  the  dissipation  in  a  collision  free 
shock  produces  a  random  distribution  of  magnetohydro¬ 
dynamic  waves.  These  waves  are  then  treated  as  the 
fundamental  particles  of  the  plasma.  A  rough  kinetic 
theory  is  developed  which  estimates  the  heat  conductior^ 
coefficient  due  to  the  waves.  Using  this  heat  conduction 
coefficient,  the  shock  thickness  is  estimated  to  be  about 
folir  times  the  characteristic  ion  Larmour  radius.  This 

A 

prediction  is  in  rough  agreement  with  experimental  re¬ 


sults  obtained  in  a  MAST  device. 


\ 
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COLLISION  FREE  MAGNETOHYDRODYNAMIC  SHOCK  WAVE* 

by 

A.  Kantrowitz,  R.  M.  Patrick  and  H.  E.  Petsch'k 

Avco -Everett  Research  Laboratory 
Everett,  Massachusetts 

Introduction 

Collisional  dissipation  in  plasma  becomes  very  slow  at 
hi(?h  temDC ratures .  It  seems  likely  that  other  dissipative  mechanisms 
associated  with  what  has  loosely  been  called  "magnetohydrodynamic 
(MHD)  turbulence"  will  then  become  of  prime  importance.  A  shock 
wave  propagating  perpendicular  to  the  magnetic  field  in  a  plasma  where 
the  ion  cyclotron  radius  is  much  smaller  than  the  mean  free  path  pro¬ 
vides  a  good  opportunity  to  study  these  phenomena.  This  paper  reports 
progress  on  theoretical  and  experimental  studies  of  these  shock  waves. 
Wave  Steepening  and  Thickest  Shock  Hypothesis 

In  many  astrophysical  or  laboratory  cases,  shock  waves 
are  formed  by  the  steepening  of  gradual  compression  fronts.  As 
the  steeoness  of  the  front  increases  diffusion  processes  become  im¬ 
portant  and  at  some  steepness  can  transfer  sufficient  momentum  and 
energy  and  produce  sufficient  entropy  so  that  a  steady  state  shock  pro¬ 
file  is  attained.  One  would  therefore  expect  that  the  diffusion  process 
which  can  act  at  the  minimum  steepness,  i.e.  at  the  longest  range,  will 
be  the  one  which  controls  the  shock  structure  provided  that  it  can  pro¬ 
duce  sufficient  dissipation.  This  expectation  will  be  referred  to  as  the 


*  Sponsored  jointly  by  the  United  States  Air  Force,  Office  of 
Scientific  Research  and  Office  of  Naval  Research  under  contracts 
AF  49(638)-6l  and  Nonr-2524(00)  respectively. 


thickest  shock  hypothesis.  + 

3  4 

It  has  been  suggested  by  Kahn  and  by  Parker  that  a 

dissipative  mechanism  can  be  found  in  terms  of  plasma  oscillations 

which  leads  to  a  shock  thickness  of  the  order  of  the  .Debye  length. 

5 

Gardner  et  al  have  suggested  that  a  permanent  shock  structure  can 

be  found  from  a  sevies  of  pulses  whose  dimensions  are  of  the  order  of 

the  characteristic  electron  gyro-radius^^^^.  where  m  is  the 

electron  mass;  N  is  the  electron  density;  e  is  the  electron  charge;  and 

c  the  velocity  of  light.  In  this  paper  we  will  attempt  to  show  that 

in  the  presence  of  a  transverse  magnetic  field  ,  randomized  MHD 

waves  can  provide  the  dissipation  necessary  for  a  steady  state  shock 

structure  based  on  a  length  fjj  the  characteristic 

ion  gyro-radius  (M  is  the  ion  mass).  The  dissipative  mechanism 

discussed  here  can  provide  the  required  dissipation  for  shocks  whose 

velocity  is  less  than  about  five  times  the  upstream  Alfven  velocity 

5).  For  non-  relativistic  plasmas  this  length  is  orders  cf 

magnitude  larger  than  the  Debye  length  or  the  electron  gyro  radius 
+  Several  examples  of  the  operation  of  the  thickest  shock  hypothesis 

in  collisional  shock  waves  can  be  cited.  For  weak  shock  waves  in 

CC>2  the  dissipative  mechanism  is  provided  by  the  lag  of  vibrational 

energy  behind  the  changes  in  translational  energy.  Steady  shock  structures 

have  been  observed  experimentally  which  are  of  the  order  of  10^  times 

thicker  than  the  mean  free  path.  ^  Similarly,  in  magnetohydrodynamic 

2 

collisional  shocks,  Marshall  has  shown  that  the  Joule  dissipation 
can  provide  the  necessary  transport  and  establish  a  thick  shock  wave 
unless  the  shock  strength  is  too  large. 
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and  thus  13  preferred  according  to  the  thickest  shock  hypothesis. 
Proposed  Collision  Free  Shock  Model 


The  model  we  propose  is  based  on  three  facts: 

1.  MHD  modes  exist  whose  group  velocity  is  faster 
than  the  propagation  velocity  of  shock  waves  (M^  5).  Thus,  these 

waves  can  transfer  energy  upstream  inside  the  shock  front.  Con¬ 
sidering  small  amplitude  waves  propagating  in  a  homogeneous. 

Cold.  CO  Hi  S  ion  -  f  r  nh«ma  am  H  n#»gl  parting  rti  «pl  a  r  ^m#»nt  pijrrpnt^ 

the  dispersion  relation  is 


(1) 


where  is  the  frequency  and  A  the  wave  number  of  the  waves, 

fj 'x  is  the  uniform  applied  magnetic  field,  and  0  is  the  angle 

between  and  k,  ,  From  eq.  (1)  the  energy  flow  perpendicular 

to  the  magnetic  field  (i.e.  the  x  component  of  the  group  velocity  V  ) 

can  be  found.  We  shall  be  interested  in  propagation  directions  chosen 

to  maximize  V  for  a  given  k.  The  maximized  V  valuesfwhich 
gx  gx  X 

generally  result  for  k  Ot  45°  to  the  magnetic  fieldfare  shown  in 
Fig.  1.  It  will  be  seen  that  energy  can  be  transferred  upstream  at 
speeds  up  to  about  five  times  the  Alfven  velocity  (  5Va>- 

2.  MHD  waves  superposed  on  a  compression  front  will 
grow  in  energy  by  compression.  A  wave  packet  will  exert  a  pressure 
of  the  order  of  magnitude  of  the  wave  energy  per  unit  volume  (<£) 
on  the  surrounding  medium.  If  the  surrounding  medium  is  undergoing 
a  compression,  the  work  done  against  the  pressure  exerted  by  the  wave 
packet  will  appear  as  increased  energy  of  the  wave  packet.  In  the  case 


— 


, 
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of  a  one -dimensional  compression  front  propagating  in  the  x  direction, 
this  work  is  given  by  “fag  per  unit  time  where  is  the  stress  exerted 

by  the  wave  packet  in  the  x  direction  across  a  plane  perpendicular  to  x. 

It  can  be  seen  that  if  a  wave  packet  spends  a  time  in  the  shock  front  equal 
to  the  particle  time  in  the  shock  front,  then  the  ratio  of  its  energy  across 
the  shock  is  similar  to  that  experienced  bv  a  c?as  undereoinp  a  similar 
isentropic  compression.  If,  however,  a  wave  packet  spends  a  long  time 
in  the  shock  front,  its  energy  can  grow  indefinitely. 

3.  The  scattering  of  waves  on  waves  can  provide  the  random¬ 
izing  process  which  is  necessary  to  increase  entropy  in  the  shock  front. 
Waves  capable  of  affecting  the  dynamics  of  a  shock  frmt  will  have 
sufficient  amplitude  so  that  non-linear  effects  will  be  important.  Thus, 
in  particular,  appreciable  changes  in  magnetic  field  strength  and  propa¬ 
gation  velocity  will  occur.  (However,  in  the  interesting  k  range,  kr^  f\J  10, 
the  accompanying  density  fluctuations  and  the  plasma  kinetic  energies  are 
small.  )  In  analogy  to  the  kinetic  theory  of  gases  we  can  estimate  a  mean 
free  path  \w  travelled  by  a  wave  packet  before  appreciable  scattering 
is  produced  by  the  random  field.  A  wave  will  be  appreciably  altered  in 
amplitude  and  direction  when  its  phase  has  been  changed  by  unity  (or  its 
"optical  path"  changed  by  1 /k)  due  to  the  disturbing  waves.  If  we  approx¬ 
imate  Eq.  (1)  by  neglecting  the  angular  dependence  of  the  dispersion 
relation  and  assuming  I  ,  it  may  be  written  as 


< a 


e  Hr.  fak)3" 


(2) 
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I 


At  constant  frequency  the  change  in  wave  number  is  then  related 
to  the  perturbing  magnetic  field  AH  by 

T  '  T  \lXb)  «■  «('+aT^/. 


The  perturbing  field  of  the  other  waves  will  be  coherent  over  a  distance 
of  the  order  of  the  reciprocal  of  the  mean  wave  number,  i.e.  . 


in  gome  this  <11  stance  tne  phase  ol  trie  wave  therefore  ^uutigsz  by  ar.; 

|f 

ammint** _ m  *r3  Th#»  nhasp  rhanffps  hv  a  rsndnm  walk  nrocess 


The  phase  changes  by  a  random  walk  process 


amount* 


with  steps  of  this  amount  each  time  the  wave  goes  a  distance  l/km. 

The  length  ^  required  to  obtain  an  r.  m.  s.  phase  change  of 

unity  (the  ne an  free  path}  is  therefore  given  by  the  number  /KK  °< 


steps  required. 


v  =  _ 1 _  Atf  <4> 

^  Ik) 

where  (ffi/ *s  ra**°  °f  average  wave  energy  to  the  field 

energy.  In  later  developments  we  will  set  =  k  thus  assuming  that 

only  waves  closely  grouped  about  a  single  value  of  k  have  appreciable 


energy. 


We  are  now  prepared  to  construct  a  rough  model  of  a  shock 


front.  Fig.  2  is  drawn  for  a  shock  propagating  at  twice  the  Alfven 
velocity  into  a  cold  plasma.  Our  present  state  of  knowledge  does  not 
allow  us  to  say  anything  about  the  shock  profile.  For  purposes  of 
illustration  and  estimation  of  the  shock  thickness  it  has  been  assumed 
that  the  flow  velocity  and  the  average  magnetic  field  vary  linearly 
with  distance  through  the  shock  front.  Since  waves  are  the  dissipative 
mechanism  the  randomized  energy  will  appear  as  wave  energy  behind 


■ 


the  shock,  i.  e. 


j£ 

fT<V 


(5) 


Variations  of  the  expression  on  the  left  evaluated  inside  the  shock  front 
will  depend  on  the  precise  nature  of  dissipative  energy  and  momentum 
transfeis  inside  the  shock  front.  At  the  present  writing  we  cannot  say 
much  about  how  this  quantity  will  vary  through  the  shock  front.  We 
will  assume  by  analogy  with  ordinary  collisional  shocks  that  it  is 
constant  through  the  shock  front.  With  the  further  assumption  that 
the  pressure  and  energy  density  cf  the  waves  are  equal,  the  conserva¬ 
tion  laws  now  permit  us  to  calculate  a  required  energy  transport 
upstream,  q,  which  is  also  plotted  non-dimensionally  in  Fig.  2.  Note 
that  in  the  center  of  a  shock  with  =  2  the  wave  energy  has  an  average 
upstream  group  velocity  about  equal  to  the  undisturbed  Alfven  speed. 

E s timate  of  the  Shock  Tnickness 

The  physical  picture  given  above  will  bf,  used  to  obtain 
order  of  magnitude  estimates  of  the  dominant  wave  length  generated 
in  the  shock  front  and  of  the  shock  thickness.  The  energy  flow  upstream 
will  be  a  diffusion  process  from  high  to  low  wave  energy  regions. 

Roughly  one  can  imagine  the  scattering  of  wave  packets  to 
be  similar  to  the  collision  of  gas  molecules.  The  gradient  of  the  wave 
energy  will  establish  an  anisotropic  distribution  of  these  waves  which 

dx 

will  result  in  a  net  energy  flow  upstream.  On  the  other  hand,  the 
scattering  processes  will  tend  to  create  sotropy. 

In  analogy  to  the  kinetic  theory  of  heat  conduction  we  will 
assume  that  the  upstream  energy  flow  is 

i  =  X*'(v”r-u5 
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(*>) 


! 


The  thickest  shock  hypothesis  can  now  be  employed  to  estimate  the 

dominant  wave  number.  Clearly  for  a  given  energy  density  gradient 

waves  with  a  maximum  /US.-*)  will  have  the  greatest  diffusivity 

and  thus  will  lead  to  the  thickest  shock  wave.  We  consider  the  case 

=  2  and  evaluate  all  quantities  (except  *1’ ♦  ^  )  at  the  center  of  the 

shock.  From  eq.  (4)  and  the  calculations  of  V  quoted  above,  the  value 

S ** 

of  k  for  maximum  diffusivity  was  found  to  be  and  the  diffusivity 

at  this  wave  number  was  *-Wa  •  The  conservation  laws  give  the 
energy  gradients  Z  .  Thus,  equating  the  energy  flow  re¬ 

quired  by  the  conservation  law  with  that  given  eq.  (6),  yields  a  shock 
thickness  L.  =Hr£. 

Similar  calculations  for  other  Mach  numbers  yield 


for  M  =  |.fT 

M  m3 


l  =  10  r 


L  =  i 


r. 


It  is  hardly  necessary  to  point  out  that  these  thickness  cal¬ 
culations  are  only  order  of  magnitude  estimates.  It  is  very  easy  to  make 
equally  plausible  arguments  for  shock  thickness  varying  a  factor  of  3 
either  way  from  the  values  given  here. 

Comparison  with  Shock  Tube  Experiments 

Shocks  propagating  perpendicular  to  an  initial  magnetic  field 

in  hydrogen  have  been  produced  at  this  laboratory  in  an  annular  electro- 

6 

magnetically  driven  shock  tube  (MAST).  With  MAST  it  is  possible  to 
calculate  the  strength  of  the  driving  field  and  thus  from  the  conservation 
laws  the  expected  shock  velocity.  The  velocity  of  the  shock  thus  calcu¬ 


lated  agrees  very  well  with  experimental  measurements.  Shock  velocities 
up  to  about  56  centimeters  per  microsecond  have  been  produced.  The 
density  ratio  across  the  shock  can  also  be  calculated  from  the  conservation 


laws.  The  density  behind  the  shock  front  can  be  found  from  brcm- 


strahlung  measurements  in  the  visible  region  with  a  photomultiplier. 
Interpreting  this  bremstrahlung  in  the  terms  of  the  plasma  density 
assumes  (1)  that  the  gas  is  fully  ionized,  (2)  chat  none  of  the  hot 
plasma  is  lost  to  the  waU,  (3)  that  no  appreciable  contamination  is 


present.  Making  these  assumptions,  the  denisties  have  been  calcu¬ 
lated  from  the  visible  radiation  and  good  agreement  with  the  conserva¬ 


tion  law  values  for  the  cases  of  magnetic  field  parallel  r  H )  and 

perpendicular  H  -  to  fhe  shock  propagation  direction  have  been 


found. 


Shocks  have  been  produced  in  which  the  characteristic  ion 
cyclotron  radius  is  as  much  as  100  times  smaller  than  the  mean  free 
path  downstream  of  the  shock  wave.  It  has  been  found  that  the  density 
(light  intensity)  jumps  and  thus  presumably  the  shock  waves  are  much 
thinner  than  the  downstream  mean  free  path.  It  is  not  clear  that  the 
light  intensity  rise  is  indeed  the  shock  thickness  since  it  might  be 
expected  that  the  shock  is  not  accurately  plane.  The  shock  thickness 
measurements  presented  on  F'g.  3  are  thus  rough  upperbounds  and 
as  can  be  seen  a  considerable  scatter  has  been  obtained. 


a  ne  order  of  magnitude  shock  thickness  estimate  previously 
made  has  also  been  plotted  on  t  ig.  3  and  is  in  fortuitously  good  agree¬ 
ment  with  the  measured  shock  thickness.  Much  work  remains  to  be  done 


before  this  agieement  can  be  taken  seriously. 
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WAVE  NUMBER  kr. 


The  x  component  of  the  group  velocity  Vgx  for  ’’fast"  MHD 
waves  propagating  in  a  cold  plasma.  The  fast  mode  is  cir^ 
cularly  polarized  in  a  direction  opposite  to  the  direction  of 
rotation  of  the  ions  in  the  undisturbed  magnetic  field  Hz, 
This  mode  is  not  heavily  damped  until  frequencies  close  to 
the  electron  cyclotron  frequency  are  reached.  The  plot  is 
made  for  waves  having  propagation  k  chosen  to  maximize 
the  VgX  for  a  given  |  k|  . 


FLUCTUATING 
FIELD  LINES  ' 


W' 


nCwWJM  l  u  * 


SHOCK 

CENTER 


^  ENERGY 
DIFFUSION 


UPSTREAM- 


i?ig. 


Physical  Picture  of  a  Collision  Free  Shock  Wave.  The 
directed  kinetic  energy  is  converted  in  the  shock  front  to 
randomized  wave  energy.  The  random  wave  amplitude 
required  is  estimated  from  the  conservation  laws.  The 
consci %'at>on  laws  and  the  steady  state  condition  also 
permit  an  estimate  of  the  average  upstream  flux  of 
energy  and  momentum  which  is  also  plotted. 
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COLLISION  SHOCK 
THICKNESS 


3  Comparison  of  theoretical  and  experimental  collision 

free  shock  thicknesses.  Thickness  measurements  from 
bremstrahlung  and  theoretical  estimates  from  the  wave 
interaction  model. 
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